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At temperatures below the magnetic anisotropy energy, monodomain magnetic systems (small 
particles, nanomagnetic devices, etc.) must relax quantum mechanically - thermal activation is 
ineffective. The discrete nature of the spectrum is important. This quantum relaxation must be 
\ mediated by the coupling to both nuclear spins and phonons (and electrons if either particle or 

. substrate is conducting). 

We analyze the effect of each of these couplings, and then combine them. Conducting systems 
can be modelled by a "giant Kondo" Hamiltonian, with nuclear spins added in as well. At low 
temperatures, even microscopic particles on a conducting substrate (containing only 10 — 50 spins) 
" will have their magnetisation frozen over millenia by a combination of electronic dissipation and 

, the "degeneracy blocking" caused by nuclear spins. Raising the temperature leads to a sudden 

' unblocking of the spin dynamics at a well defined temperature. 

Insulating systems are quite different. The relaxation is strongly enhanced by the coupling to 
nuclear spins. At short times the magnetisation of an ensemble of particles relaxes logarithmically 
^ ' in time, after an initial very fast decay - this relaxation proceeds entirely via the nuclear spins. At 

longer times phonons take over, but the decay rate is still governed by the temperature-dependent 
nuclear bias field acting on the particles - decay may be exponential or power-law depending on the 
temperature. 

The most surprising feature of the results is the pivotal role played by the nuclear spins. The re- 
1^-^ ' suits are relevant to any experiments on magnetic particles in which interparticle dipolar interactions 

I are unimportant. They are also relevant to future magnetic device technology. 

^ '• PACS numbers: 75.10.Jm, 75.60.Jp 

s: 
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A. Quantum Relaxation 



One of the most thoroughly explored subjects in all of science is that of thermal relaxation of magnetisation, in 
magnetic systems of all shapes and sizes. Investigations in this area (which go back many centuries |^) have revealed 
^ I many subtleties, and even today there are many unsolved puzzles (e.g., the physics of "magnetic avalanches"). 
- - ' Very recently a whole new set of questions in this area has arisen, with the advent of well-characterised " nanoscopic" 
magnetic structures Such structures include " made-to-order" magnetic grains, magnetic wires and superlattice 
arrays, as well as thin films and spin chains. There are very obvious applications of such nanomagnets in, e.g., the 
recording industry, as well as in information transmission and computing. Theoretical activity has been particularly 
intense on 1-dimensional spin systems j^], and on the tunneling and coherence phenomenon which should exist in 
both grains (involving "giant spin" dynamics) and in the dynamics of domain walls Both the theory and the 
experimental activity in the latter field have been reviewed recently . 

In this paper we focus on something rather different - the " quantum relaxation" of magnetisation that must occur 
in any small magnetic particle, once thermally-activated relaxation has ceased. Although there does not seem to have 
been a serious theoretical analysis of such relaxation before, the conventional view has been that it must proceed via 
a simple tunneling process (dissipative or otherwise). The picture we shall derive here is more subtle, and contains 
some rather fascinating new physics. Tunneling is certainly involved; but we find that what really controls relaxation 
is the nuclear spin system inside the magnetic particle, often with the help of phonons or electrons. One can put the 
result in the following way: the nuclear spins provide the essential "switch" that turns on (or turns off) the quantum 
transitions. As far as we are aware, there has been no previous recognition that nuclei had any role whatsoever to play 
in magnetic relaxation, despite the enormous number of experiments in this area. Nor do we know of any theoretical 
analysis of the influence of electrons on magnetic relaxation. There has certainly been no attempt at a theory dealing 
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with the combined effects of nuclear spins, phonons, and electrons, which is what we do here. As we shall see, it is 
impossible to give a realistic analysis without putting them all together. 

We shall not try here to address experiments in any detail, for two reasons. First, each different experimental system 
has its own idiosyncracies, and we shall see that our theory has many different parameters - it is simply not possible 
to give detailed results for all possible cases, and we feel that it would be preferable to examine particular systems 
for which parameters are well known, as test cases. This we intend to do elsewhere. Second, in most experiments on 
magnetic grains, the grains are rather close together. As we shall briefly explain in this paper, this makes it probable 
that even the low-T relaxation of an experimental ensemble of grains proceeds via dipolar interactions between the 
grains. 

The theory in this paper is closely related to previous work on quantum coherence in magnetic grains [ p^ . In that 
work a very detailed analysis was given of the way nuclear spins control the dynamics of "giant spins", i.e., spins 
with quantum number S ^ 1. We also analysed the effect of phonons and electrons in "decohering" the quantum 
motion of S. The problem of quantum relaxation is different, in that even if coherence is destroyed, relaxation may 
still proceed incoherently. What we shall find is that in the absence of nuclear spins, phonons cannot relax S* in a 
small biasing field, because they couple too weakly to it; on the other hand electrons couple so strongly to S that they 
freeze it completely. The nuclear spins can then liberate S - amusingly, the mechanisms by which they do this are the 
same as those by which they destroy the quantum coherent motion of S. The bias field e acting on S, generated by 
the combined hyperfine fields of all the nuclei can easily be greater than lOi^T (e.g., for TbFes grains), and depending 
on the size of the grain may be equivalent to an external field from 100 G up to well over 1 Tl This bias changes 
with time in a diffusive way, at a rate governed by inter-nuclear dipolar interactions (in the absence of phonons and 
electrons), and also nuclear spin-phonon or spin-electron interactions, when the latter are present. Every so often, as 
e shunts around, it brings quantum 2 states corresponding to different orientations of S into near degeneracy. At this 
point the system can tunnel, provided (a) the nuclear bias field e(t) gives it time to do so coherently, (b) the usual 
dissipative effects of the electrons and phonons are not too destructive. 

On the other hand, we shall also see that the nuclear spins can also further inhibit the relaxation of S, at low T. If 
the nuclear spins are allowed to come to equilibrium, then at temperature below the hyperfine coupling energy, they 
exert a large negative bias field on S, which essentially traps S for astronomical times. Thus, depending on T and on 
the nuclear polarisation, the nuclei can enhance or suppress relaxation. 

Although our purpose in this paper is not to deal with the implications for magnetic technology involving quantum 
relaxation, we feel we should at least make some remarks here about what they might be. Fig.l shows the way in 
which magnetic computer memory elements have decreased exponentially in size in the last 40 years - such plots have 
figured large in recent discussions of the need for new kinds of computer design, incorporating elements which operate 
quantum-mechanically. At present computers, as well as magnetic tapes, use elements which behave classically, and 
are stable over long periods of time. This stability over decades exists because they are big - the energy barrier 
between 2 states of the element, proportional to its size, usually exceeds 100 fc^T even at room temperature. The 
hypothesized "quantum threshold", below which tunneling is important even at T = 0, is usually supposed to be for 
grains containing roughly 10^ — 10'* spins, depending on the material involved. However we shall find that the threshold 
between stable behaviour and quantum relaxation can be moved around a lot, depending on how the magnetic system 
is coupled to the nuclear spins, and whether or not it is coupled to electrons. We shall see that in some cases, one can 
freeze the dynamics of even microscopic magnetic systems (containing only 10 spins), for times of thousands of years 
or more. Conversely, we shall see how it is possible for much larger systems (containing 10^ electronic spins or more) 
to relax very quickly (in /is). These theoretical insights may have some relevance to the design of future nanoscopic 
magnetic devices. 

In what follows we shall analyse the role of electrons (Section ||), phonons (Section pl| ), and nuclei (Section |rv| ), 
acting in isolation on S. Then in Section we shall put them together, to give the final rather complex picture. 
As noted above, we save detailed discussion of experiments for other papers; but in Section VI we will indicate the 
general nature of our predictions. 



B. The Model, and Energy Scales 

We briefly describe the giant spin model here, and the various energy scales relevant to the physics. The basic idea 
behind this model has been explained in detail elsewhere |9|JTl| . The exchange couplings Jy between electronic spins 
in the grain are enormous compared to the anisotropy couplings Ka (typically Jij eV, whereas Ka ~ 10^^ — 1 K). 
The mesoscopic grains can be treated as giant rotators, keeping only states | S, m), with S > m > —S. The spectrum 
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of this giant spin is shown in Fig.l in a small bias field ^h- Sometimes, as a model example, we shall use the 
easy-axis/easy-plane system, having Hamiltonian 



S 



K« Si + K^Sl 



(1.1) 



We will usually assume Ho = zHq, so S^h — jSHo- 

The low energy physics of (LI) can be understood in terms of the truncated Hamiltonian 



Ho{t) = {2Aofx costtS - ^hTz) = '^sTx ~ S.hTz , 
where the tunneling splitting Aq is 



and the bare parameters are 



Ao rJoe 



(1.2) 
(1.3) 
(1.4) 



Ao-^2SiK\\/K^y/^ , (1.5) 

The "bounce frequency" ilo is also roughly the distance to the next pair of excited levels; the total barrier height 
in (pTll), between the 2 semiclassical minima at S = ±zS, is SK\\. Thus if, say, S ~ lO'*, we might easily have 
flo 0.1 - 1 if, a barrier height ~ 0.01 - 0.1 eV, but a splitting Aq ^ 10"^ K, or less. 

This is certainly not the whole story. Magnetic grains having more complicated symmetries are not so easily 
described using WKB/instanton methods. The "internal" magnon modes of the giant spin are also neglected, except 
for those uniform precession modes corresponding to | 5, m) . There will also be phonons (spectrally very weak at 
these low energies) and possibly electrons (which will mix strongly with all states of the grain). 

Finally and most importantly, the nuclei in the grain (and also possibly outside) will couple via hyperfine interactions 
to each level shown in Fig. 2. The net result of this is shown in Fig. 3 for the 2 lowest levels in the presence of a bias 
field. Typically f2o ^ ^ Aq, where cok is an individual hyperfine coupling; however N nuclei will spread out 
a single grain level into a mass of 2^ levels, of Gaussian half-width ~ N^^'^ujo, where lUo is the principal hyperfine 
splitting. For rare earths, where Wo ~ 1 GHz or greater, this half-width can easily exceed Qq for mesoscopic grains, 
and the whole "giant spin" model begins to get rather complex. 

In this paper we will ignore all physics at energies ~ f2o or greater, simply including it into renormalisation 
parameters of a low-T effective Hamiltonian. For such a treatment to be valid, we then require 

ksT, (h, N^'^uJo « ^o . (1.6) 

We will also assume that any "loose spins" on the surface of the grain still have couplings ^ fio to S. If not, they 
will be treated as part of the "spin bath" environment, along with nuclei and paramagnetic impurity spins. 

Before beginning, it is useful to bear in mind 2 simple points. Note first that in the presence of any bias S^h which 



is greater than Ag in (1.2), the system is essentially "frozen" in one or other of the states Iff) or |-|J.). Suppose, e.g., 
the system starts in the higher-energy state Iff) at t = 0. Then the probability p(°^(t,^//) of finding it Ifl") at time i, 
in the absence of any other couplings, is given by 

p(°)(t,eH) = l-|f sin^iJi, (1.7) 

where E'^ = ^ ^s- Thus the coupling to phonons, electrons, and nuclei is necessary if we are to have any 
relaxation at all. Usually one thinks of this as a kind of spontaneous or stimulated emission, perhaps in the presence 
of dissipation, in which the requisite energy is taken up by a bath of oscillators. 

In fact this conventional picture of tunneling of a biased 2-level system, coupled dissipatively to a bath of oscillators 
[ p^ , works reasonably well for a giant spin coupled to either electrons or phonons. This we shall see in detail in the 
next two sections. However as soon as we introduce the spin bath of nuclei, etc., it breakes down completely; this is 
basically why our final results are so surprising. 

The second point arises from the finite-dimensional nature of the giant spin Hilbert space. As emphasized by van 
Hemmen and Suto [Hjl^ this means that we must use WKB calculations of tunneling rates with great care when 
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discussing the motion of S. This point is especially clear in recent work of Politi et al. |14| (see also Villain et al. 
0), analysing the Mni20i2 acetate systems; here tunneling appears to proceed via a 4th-order term ~ (S^ + St) 
in Ho{S). Consequently transitions may only proceed between states differing by ASz — ±4. To correctly handle 
such selection rules in the instanton calculus is not a trivial problem. Even more important, there is no continuum 
of final states in this tunneling problem, and tunneling may only proceed if there is near resonance between initial 
and final states. This means that in the presence of a bath of, e.g., phonons, one may not in general simply apply the 
Caldeira-Leggett tunneling formalism [ p5[ to calculate tunneling, since this formalism assumes a continuum of final 
states. For this reason calculations such as those of Garg and Kim of the tunneling rate of a strongly biased grain 
interacting with phonons, can only be correct when S oo. For finite S one must take account of the discrete nature 
of the spectrum of S, at least when the level spacing is greater than the (very small) linewidth caused by phonon 
damping (note further that calculations of grain tunneling in the presence of a nuclear spin bath should never use the 
Caldeira-Leggett method, since the spin bath cannot be mapped onto an oscillator bath of the Feynman- Vernon jl^] 
type - thus calculations of this kind are not valid for any S). 

We now turn to the detailed treatment of the different relaxation mechanisms. Readers wishing to see only a 



summary of the results should proceed directly to section VI 



II. GIANT SPINS AND THE ELECTRON BATH 



In this section we develop a model for the coupling between the giant spin and the electron bath only, and use it 
to calculate relaxation rates for the grain magnetisation. Since, as far as we are aware, no attempt has ever been 
made to deal with the coupled giant spin/electron system before, it is clear that we must begin from first principles. 
The full development of our model, including detailed discussion of the multiple-scattering of the electrons by the 
individual spins inside the grain, the effects due to surface spins, to the discreteness of the electronic spectrum inside 
the grain if it is on an insulating substrate, the interaction between electrons and the internal degrees of freedom 
(magnons) in the grain, etc., will not appear in this paper. We shall also set aside questions about how conduction 
electron dynamics outside the grain, in a conducting substrate, will be influenced by the giant spin dynamics. This is 
a very important question for any experiment (which can use the substrate magnetoresistance as a probe of the grain 
dynamics), but is too complicated to deal with here. 

There are 3 physical situations which we will address here. The first involves conducting electrons both inside 
and outside the grain; we shall see that this problem can be analysed fairly clearly. Likewise the problem of an 
insulating grain interacting with conducting substrate can be modelled in a reasonably treatable way. The problem 
of a conducting grain and an insulating substrate is more messy, however we can fairly easily extract the essential 
results. Our essential result is that unless the substrate is insulating, electrons simply freeze S , unless S is very small. 
We emphasize that in reality, nuclear spins must be included (Sections IV and M). 



A. Conducting Grain and Conducting Substrate. 

We shall start by assuming that the giant spin, contained within a volume Vo ^ R^, interacts with an electron fluid 
which permeates freely through the boundary between grain and substrate; i.e., a Hamiltonian 

H = H.is) + c^^^c,^^ + ^ E • E E ^'H^.^A ' (2-1) 

fe,cr le^o k q 

describing short-range exchange interaction between the localised individual spins si, located at positions ri inside the 
grain, and conduction electrons in momentum \ k). A more refined model would couple grain electrons in eigenstates 
I /z) to substrate electrons | k) , via a transfer matrix ^ across the boundary. However it is intuitively obvious that 
this is unnecessary once | T^j^ |> Ae^, where Ae^ is the typical level spacing between internal electron states of the 
grain. If the electronic bandwidth is D, then Ae^ ^ D/S, and for most metals D ^ 1 — 5 eV; thus if S* ~ 10^, 
Ae^ 10 — 100 K, so that an insulating layer at least 10 A thick would be required to make such a model necessary. 

We assume in using ( |2.l| ) that the internal grain spins are locked together to form a giant spin, and that the Ji 
already take account of the virtual mixing between internal excited states of the magnetic ions and of the giant spin 
(magnons, surface modes) and the electron bath. Thus, as usual, we shall be working at low T and small Ho- However 
we shall assume that the Ji can vary; near the surface of the grain boundary they may be a little weaker. The present 
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model does not include "loose spins", i.e., surface spins which are more weakly coupled to the motion of S (the reader 
should not confuse the Kondo coupling Ji here with the inter-spin exchange Jij between spins inside the grain) . These 
will be discussed in section [V onwards. Typically the Kondo couplings Ji '--^ 0.1 — 1 eV. 

We now rewrite (pj), incorporating the giant spin hypothesis, to get a volume averaged interaction 



Hi, 



/-^ Z-^ 1 k+q,a k4 



(2.2) 



/ (2-3) 

where the form factor integrates the number density p(f) of microscopic spins over the grain volume, and J is the 
mean value of the Ji. is a sort of "giant spin Kondo Hamiltonian" , and as such contains a lot of interesting 

physics, much of which we will ignore in this paper. Recall that the Kondo coupling for microscopic spins is just 



( p. 21) with p{r) in (2^) restricted to a single lattice cell; for this problem one conventionally defines a dimensionless 
coupling g = JN{0) ~ 0.1 for most metals, where iV(0) ~ p/D is the Fermi surface density of states. The Kondo 
problem can be mapped onto an Ohmic spin-boson problem of the kind discussed by Leggett et al. |l2j ; even for this 
spin-1/2 system, there is a large dimensionless Ohmic coupling to the electron bath; uk ~ 0{\). 



It is possible to start with (2.1) and build up an analogous description by considering the individual Kondo scattering 
from each spin, in many different orbital channels. Here we shall finesse entirely this move by instead integrating out 
the electrons directly, starting from (^|^). We first truncate directly to a biased spin-boson Hamiltonian of form (for 
Ho{S) = H^{S) in (P)): 



^ N N 
^eff = ^oTx - + - ^ mk{xl + UjIxI) +Tz^ CfeXfe , (2.4) 
fc=l k=l 

where the bias field is parallel to the easy axis (we shall ignore transverse fields here). The Caldeira-Leggett spectral 
function for this problem is 

and the Ohmic coupling constant af, is determined by a standard Fermi surface average [ p^|j2C| ] over electrons, i.e., 

a, = 2g^S^\F^^^j:,\'')F,s.. (2.6) 
The naive guess that af, ~ 5^, apparently confirmed by (^.(j), is however wrong; carrying out the average, one gets 

This at first surprising result is a consequence of interference between Kondo scattering at different sites - notice 
that k'^^ <C i?o, so that the problem is a multiple scattering one. A related point is that the estimate of a.h in 
( ^.7[ ) looks perturbative - it is not clear it should work for large a^. However for a grain of radius Ro some Z^n^^. ^ 
l/{kFRoY ^ 5'^/'^ orbital channels are playing role in electron scattering. If S" < 10'^, the phase shifts in each channel 
are small (one has 5i ^ gkpRo ~ 9^^^^), making it possible to apply a perturbation theory, yet the sum over all 
the channels can be rather large. The case of large phase shifts is more complicated, and may result in a giant spin 
dynamics different from that described by the standard Ohmic model. We address this possibility in Appendix A. 

For mesoscopic spins {S > 10'^) the coupling in (2.7) is very large; we may take over standard results from the 



spin-boson problem to discuss its effect on the dynamics of S. Since a?, > 1, coherence at zero bias is destroyed at any 
temperature; at T = the grain dynamics are rigidly locked, and the spin S cannot move. At finite T and e = 0, the 
system fluctuates incoherently between the states Iff) and at a rate ~ (r/Oo)^"'"^; on the other hand if e ^ T, 
one has incoherent relaxation at a rate ~ (e/Jlo)^"''"^. The general result for the relaxation rate is p^,pO[ 
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^"^^ / ^ \ I I irVL -1- If: I 'ATT /II- 

(2.8) 



where r[z] is the Gamma fmiction; for e = 0, this gives 



, ) r[2a,,] ) 
On the other hand for very low T, such that e/kBT ^ 1, it gives 



-iT).^A^Y^(^-^T-\ (2.9) 



re-^W = 2vrAj'^)-i-^f-lV""- (2.10) 



no) r[2ab] 

The crucial thing to notice here is that for microscopic spins (5 ^ 100), there will be easily observable relaxation of 
5. On the other hand if S is much greater than 100, the giant spin will be frozen completely, even over astronomical 
times, unless either bias or temperature is large (e, ksT ^ fio). Moreover, in this model we would expect that applying 
either a large bias or temperature would lead to a sudden "switching on" of the relaxation at some critical bias or 
temperature (however recall that our model is not strictly valid for e, ksT ^ fi^, since higher levels of S will then 
enter into the picture). 

To see this, consider 2 examples. First, suppose that ~ 10^^ and S ~ 30, so af, w 1. Then ( p.9D shows that for 
bias e = 0, we have r-i(r) - A^T/fJ^^ and if ksT = 0, we have t-^{T) - A^e/fJ^. Since A^/fJ^ might be something 
like 100 Hz, this implies low T or low e relaxation times of fractions of a second. 

On the other hand, consider a second example where ~ 10~^ but S ^ 300, so ab ~ 20. Then if e = 0, we find 
that r-i(r) ~ 10-12 /^l/Vto {2TrkBT/no)^°, and if T = 0, we find T^^e) - lO'^e Al/no (e/f]o)'*o.Then if e = 
we see that the relaxation suddenly switches on at a temperature ksT ~ r2o/27r (if we are interested in a relaxation 
time of order seconds); for temperatures below this, relaxation times increase catastrophically (if T ~ Tc/2, we find 
T~^{T) ^ 10^ years!). The same happens as a function of bias when T ^ e. Obviously if we go to even larger grains, 
with S > 1000, this switching process becomes almost a step function of temperature. Although our 2-level model 
starts to break down once irksT, e ^ fto, this will not affect the validity of this switching result, or our calculation of 
the switching temperature. 

We have not included nuclear spins in these calculations - in section ^ we shall see how a combination of nuclear 
spins and electrons affects the dynamics of S. 



B. Insulating Grain and Conducting Substrate. 

The insulating grain can be analysed in a similar way, provided we gloss over some of the complications arising 
from the variation of both the electron density and the spin properties in the vicinity of the surface. Thus we naively 
extend the previous model by writing an effective interaction Hamiltonian 

HL^ - • E E ; (2-11) 



/ 1? I*e(r)pp(r>''", (2.12) 

where now the form factor takes account of the decay of the electron density | ^e('^) P as one penetrates into the 
sample. By the same manoeuvres as before this leads to an Ohmic coupling 

as g'sS'/' , (2.13) 

where gs ^ g, but is multiplied by some "scale factor" which describes the depth to which the electrons penetrate 
into the grain; clearly it can vary widely. 



The reduction to only near surface spins makes as <C at for large S, but from (2.13) we see that as > 0(1) if 
S > 10'^. For smaller grains, i.e., S < 100, one can have as < 1, and then results different from the "frozen spin" 



behaviour of (2.S)-(2.1C) are obtained. We may define a renormalised splitting As, and a damping rate Fs, by 
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(2.14) 



Ts = 2TrasT . (2.15) 

At T = one can get [Q coherent motion if as < 1/2; however any such coherence is destroyed once Tg > A^, 
which even for microscopic values of S (i.e., S ~ 10) wiU occur at extremely low temperatures, usually much less then 
1 mK. For higher T we find, when e = and as < 1, that r^^ is still given by ( ^.9| ), which in terms of A^ is 



T-^(T) = 2As— ^ — . (2.16) 



For finite e and T, the form (2.8) is still valid, as is (2.1C) when e ^ T. 

From all of these results we conclude that for low T and e, the relaxation dynamics of mesoscopic grains will be 
extremely slow in the presence of a conducting substrate. Thus, although our models have certainly not included all 
possible nuances of the surface physics, one conclusion is very clear and unlikely to be affected by further refinements, 
i.e., that for T, e <C J7o, mesoscopic magnetic grains on conducting substrates will have their spin dynamics completely 
frozen. However once either T or e is of order Qq or greater, these dynamics will rapidly be liberated. Again, these 
statements will be modified if nuclear spin effects are important (Section ^). 



C. Conducting Grain and Insulating Substrate. 

In the opposite extreme when the substrate is insulating, the character of the electron states is quite different - they 
form a set of discrete states | /i,(7) inside the grain, with mean spacing Ae^ ~ D/S for bandwidth D. We describe 
these states via creation and annihilation operators Aj^^ ^,A^,a-; these states will be rather complicated, even in the 
absence of spin-orbit coupling. It is very useful to separate the exchange scattering interaction into bulk and surface 
parts, where now the bulk part sums over the whole sample, assuming that all the Ji are the same, whilst the surface 
part describes the difference between H^^l^ and the true interaction Hamiltonian Hint\ thus if 

H^nt ^\Y. Mk ■ ^'^^)cUc,,f, , (2.17) 

with the cl ^ creating electron on site i, we write 
Tjbuik 1^ i^/^N t 

= ^"') E ic.c,, - Jis'- ^"')E4,«^.^/3 > (2-18) 

HL = ^ E (-^i - ^)(*^' ■ '^"'')4aCi,/3 . (2.19) 

We notice that the bulk term is diagonal in both the site representation and the exact eigenstate representation, 
because of completeness of states. This is important, because the condition J ^ flo means that the electron spins 



and the surface term 



inside the system will rotate rigidly with S, and from the last form in (2.18) we see that there is no way this adiabatic 
rotation can excite the low energy orbital states | /i, a) of the system. Thus the effect of H^^l^ is simply to renormalise 
the moment of inertia of the giant spin. 

On the other hand the surface term does have off-diagonal matrix elements (/i | H^^^ \ /i') and this may result 
in infrared electron-hole pair production when S fiips. The calculation of the corresponding coupling parameter 
proceeds in exactly the same way as in the previous subsection (the corrections to the plane-wave calculation are 
small in a large grain with kpRo ^ 1). Thus at high temperatures the internal conduction electrons constitute an 
Ohmic bath with 

as ~ glS^"" , (2.20) 
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(here gs is given by (Ji — J)N{Q) averaged over the surface). Now we notice that this result holds only at T > Ae^. 
At lower temperature the effective coupling constant goes to zero exponentially with T: 

a'J^ ^ase-^'^/^;, (2.21) 

and the infrared renormalisation of the tunneling rate stops at 

A, = Aj j ; , {for Ae^ « n^) . (2.22) 

If Ae^ ^ flo one has Ag « Aq, and we do not have to worry about conduction electrons in the grain at all because 
they will rotate adiabatically with S. 

The onset of coherent motion in zero bias can be found by comparing the damping rate Tg — 2TTa''JfT with A^. 
Expressing all the parameters in terms of S we find the criterion of coherence in the form 

T T 1 

For experiments at mK temperatures and small bias e < Ae^, this means that electronic dissipation effects can be 
ignored in such grains provided S < 10^. On the other hand for T and/or e > Ae^, we return to the case of strong 
Ohmic dissipation which looks very much the same as that for an insulating grain on a metallic substrate. 



III. GIANT SPINS AND THE PHONON BATH 



In the absence of electrons (and of nuclear spins) relaxation in a biased system proceeds via phonons. This is a 
very slow process - typical relaxation times are of order months. A number of previous studies have appeared; Garg 
& Kim [|6j and Chudnovsky |^ have discussed grain tunneling within a straightforward Caldeira-Leggett approach, 
assuming a continuum of final states available for the tunneling - we have already explained problems with this 
approach in the introduction. Politi et al. jl^] have given a very thorough discussion of the effects of phonons on 
tunneling in the Mni20i2 acetate system, including the effects of discrete initial and final states, and the non-diagonal 
couplings (to be defined below) . They also take proper account of the symmetries and selection rules in the tunneling 
transitions. Whether their theory explains the relaxation measurements in Mni20i2 acetate [ p2[ is another matter 
- in our opinion the explanation may have more to do with dipolar interactions between the molecules [ p3[ |. The 
treatment of phonon effects we shall give is more general although quite close to theirs. We first discuss the form the 
low-energy effective Hamiltonian should take, depending on the symmetries of Ho{S) and the direction of the applied 
bias field, and determine the relaxation of the giant spin dynamics for the different possible couplings. The values 
of the relevant couplings are then estimated - no attempt is made at exact calculations, since these would depend 
so much on the detailed nature of the system in question. Our main result is that these couplings are so weak that 
phonons by themselves are incapable of causing anything but very slow relaxation of S. 



A. Effective Hamiltonians, and Relaxation Rates 



We start by assuming, for ease of discussion, a simple easy-axis, easy-plane bare giant spin Hamiltonian in a bias 
(Eqs.(0), ^)). The phonons are described by = E?-i'9(^9A-+l/2), with = cq (we take account of polarization 
later, and ignore optical modes), and the general form of the truncated interaction is then 

H,r.t =T.Y1 Cl{h)xg-+ [r+^Y. C^^g+H.c] ; (3.1) 
q ° q 



{2mqv^yl-^ 



(3-2) 



where h = ^eSH^ is the transverse bias. 

It is unusual to consider non-diagonal couplings (i.e., couplings in f±), since in general their effects are reduced, 
compared to those of diagonal properties, by a factor (Ao/f^o)^- However we shall see that in this system the diagonal 
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couplings are zero unless we apply a field perpendicular to the easy axis. The reason for this is very simple - time 
reversal symmetry implies that the magnetoacoustic interaction cannot distinguish between grain states Si and S2 
when Si = —S2 (a point which also arises in the interaction of phonons with domain walls ||]), and diagonal couplings 
describe interactions which do distinguish such states. Applying a transverse field projects a component of Si on 52, 
and enables such a coupling; from this argument we see that C^{h) is proportional to the transverse field H^. We 
shall later see that 



Cl^i-f,H,/K\\)C^ , (3.3) 



for a field applied in the easy plane. 



Before deriving (3.1) and evaluating the couplings, we first give the results for relaxation rates that derive from it. 
Just as for the electron bath, these rates are very low, but for the opposite reason - instead of being frozen by a very 
strong coupling, S is frozen because the effective coupling is so weak! 

The relaxation rates are derived by first defining Caldeira-Leggett spectral functions 



J„(c^) = ^V^^<5(a;-^,-), (3.4) 



2 ^ 



q 

where a = 2;,_L; these describe the phonon effect on grain dynamics via 2nd order perturbation theory. It has of 
course been known for decades that these couplings are very weak for w <C 6_d , the Debye frequency; in fact 

J^{Lo)^B^{Lo/QDr , (3.5) 

with m > 3 in 3 dimensions; this is the fundamental reason for the coherent motion of polarons and defects in 
insulators at low temperatures |2Q]. We will see below that m = 3 in our case. 



We start with the case where there is no transverse field H^j., so that C5 = in (3.1). The relaxation dynamics for 



this case are described as damped motion between S = ±£5* with a relaxation rate (see, e.g., |2^) 

W -'^^i J^^^) coth(/3£;/2) ^ AB^^ ' (3.6) 



= Ai, + e\ (3.7) 

where the longitudinal bias e — jeSHz- This rate is of course very small {t~^ ^ Ao,e). Notice that T~^{e) is quite 
different from that rate which would be calculated via a naive application of the spin-boson results Q; in fact it is 
smaller by a factor {e/flo)"^- 

On the other hand application of a transverse field H = xHx gives a diagonal coupling, according to ( |3.3|) , which 
can be incorporated into a spin-boson model; in the absence of any bias e one finds a relaxation rate 



2J,(A„)coth(/3Aj2)~4B,^ , (3.8) 



or, in a non-zero bias. 



'{h,e) « 2^ ME) coth{pE/2) ^ AB.{h)^' , (3.9) 

E^ e_D Bd 



which is independent of e as long as e <C T. Using (3.3) we may express B^ in terms of Bj^, and write 

These rates are very long, because there are so few phonons available to absorb the tunneling energy at these 
frequencies - we shall give typical numbers once the coupling constant B± has been determined. 
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B. Magnetoacoustic Couplings 



The magnetoacoustic coupling is of course very well understood |l|]; its effect on the WKB tunneling action of a 
grain has been analysed by Garg & Kim |l^ and by Chudnovsky pl| . There has been some confusion about the 
relation between Chudnovsky's " angular momentum coupling" and the usual magnetoacoustic coupling, which we try 
to resolve here. 

The essential points we wish to make are that the couplings C^, in ( |3.l| ), and B± in (|]^), are given to within 
dimensionless constants ^ 0(1) by 

\C^\^2Sno\q\; (3.11) 



^ S^nl/Qo , (3.12) 
so that the relaxation rate in a longitudinal field Hz is finally 



whereas the relaxation rate in a non-zero transverse field Hx has an additional term 

By comparing ( [3. IS ) and ( |3.14 ) we find their ratio to be ~ (iJ^SA'" / Hx^of' ■ This ratio is large unless the external 
field is almost parallel to x. 

To get some ideas of the relaxation rates implied by these formulae, let us consider a situation where S = 10^, 
Ao 1 MHz, r ~ 50 mK, Qd ~ 100 K, and e ~ "feSH^ ~ 50 mK. Then we get T(e) - lO'' sec (i.e., 3 hours). Thus, 
if we ignore nuclear spin effects, the relaxation times at 50 mK are already very long; and Aq ~ 1 MHz is actually a 
large value for Aq. If the bias energy is further reduced to the Ao scale we get astronomically large r. 

To get (3.11) and ( 3.12| ) we start by noting that Chudnovsky's angular momentum coupling 

Ham = • ^ "(^^)) ' (3.15) 

integrated over the grain volume, is nothing but that set of transverse terms in the magnetoacoustic coupling which 
is responsible for transporting angular momentum to the transverse phonons. This magnetoacoustic coupling has the 
general form 

Hme = ile/S)'^ / dfaijkiSiSjUki{r) , (3.16) 
Jg 

where the coefficients j'^aijki are energy densities, and in fact are simply combinations of the magnetic anisotropy 
constants (i.e., of if" and K± in the case of H;^{S)), since this is the only energy scale that can be involved. From 
this we see that for the grain, Hme ~ SK \ q\. However we also have 

S = i[S, HoiS)] ~ rtoS , (3.17) 

and so Ham ^ Sflo \ q \ also; the direct relation to the aijki can be easily determined for any particular Ho{S) and 

lattice type. 

From this it follows that ~ 2Sfto | g | in (3J). It is instructive to rederive this by estimating the correction 6 A 
to the bare tunneling action in a given strain field Uki, i.e., finding 

Aoiiuki}) - r!oe-^«"-» ^ A„e'^^«"-» , (3.18) 

This is determined trivially from the interaction Ham by integrating over time (assuming a negligible change in the 
instanton path); writing 
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6A{{uki}) 



we get 



df 



(3.19) 



= 25r!o / — e''«-^"(<f X %) • z 
= 2SnoFif{q X eg-) • z , 



(3.20) 



where -Fg- is the form factor, and is the polarization of the relevant phonon; for long wavelength phonons {qRo ^ 1) 



this gives us equation (3.11); and B± in (3.12) immediately follows using (3.4). Notice that in the case of the 
easy- axis/easy-plane Hamiltonian H^{S), the same argument gives 



CI 



2sno 



(3.21) 



It is sometimes useful to write all these couplings in terms of the relevant sound velocity, using ~ pc^; usually 
this will mean the transverse sound velocity since it is usually lower than the longitudinal one. Then without specifying 
exactly which of the many combinations of the aijki enter into our expression for some particular symmetry, but just 
calling this combination Ka, we write 



(3.22) 



Finally, we remind the reader that none of the analysis in this section takes account of nuclear spins. 



IV. GIANT SPINS AND THE SPIN BATH 



We now turn to a very different kind of bath from the phonon and electron bath, viz., the "spin bath" made up 
from nuclear spins both inside and outside the grain, paramagnetic electronic impurity spins, and any loose spins in or 
near the grain surface. As noted in the introduction, this spin bath cannot be mapped onto an oscillator bath, and so 
a completely different kind of theoretical framework is necessary to describe its effects on mesoscopic or macroscopic 
quantum objects. 

The main point of this section is to show how the spin bath, interacting in isolation with a giant spin S, will 
basically block any relaxation at all, unless either the grain happens to be in a total field (produced by the sum of the 
external field and the "internal bias field" generated by the nuclei) which is very small, less than roughly Ao; or the 
nuclear spin diffusion mechanism of energy bias diffusion allows the system to find the resonance window. In what 
follows wc first briefly recall the 4 mechanisms by which the spin bath controls the giant spin dynamics, and set up 
the formalism required to analyse magnetic relaxation in the presence of nuclear spins. We then derive the relaxation 
dynamics of S*, first ignoring the internuclear spin diffusion (to see how S is frozen), then including it (to show how 
S can then relax). 



A. Effective Hamiltonian, and the 4 Coherence-Blocking Mechanisms. 



In an earlier series of papers [|10|,^-26| we have constructed a theoretical description of a giant spin interacting 
with a more or less arbitrary spin bath, and solved this model for the quantum dynamics of S in the absence of any 
external bias. The same sort of model has also been set up to describe an arbitrary unbiased macroscopic coordinate 
interacting with a spin bath. Elsewhere we have used this model to deal with coherence experiments in SQUID's [p7| , 
and also shown how it leads to some remarkable conclusions about the quantum theory of measurement (particularly 
concerning the role that nuclear spins play in making the mesoscopic world behave classically p8| ) . 



11 



Here we generalise this work to deal with a biased giant spin. As before, we shall assume that the giant spin can 
be truncated to 2 levels, and that T, e <C flo', sll effects of higher levels will be absorbed into the parameters of an 
effective Hamiltonian. This Hamiltonian was derived in Refs. |lO|,p4|-p6[ , and is written for the biased case as 



N ^ N N 

Heff = 2Ao{f+ COS [$ + ^{akfik - i^kVk) ■ (fk] + H.c.} - -feS ■ Ho + Y^^k^k ■ (fk + irik ■ ■ (4.1) 

fc=i fe=i fc=i 



There is also an inter- nuclear interaction of dipolar form; thus in general we have to add to (4.1) a term 



Vi{<Jk}) - Kk^'^k^t ■ (4-2) 

We shall not need to write down the matrix elements V^j.^ in detail; their strength is ^ T2^, the transverse nuclear 
relaxation rate. Typically we shall assume that Ho is parallel to the easy axis of our easy-axis/easy-plane system, so 
that the total effect of Hg will be written as 

J J ■ Ho ^HT, , (4.3) 

where £,h = JeSHo (we change notation - the symbol "e" will be used to denote the "internal" bias coming from 
the nuclei). The spin bath is described in ( [4.1| ) by spin-1/2 variables {ct^}. The vectors Hk, Vk, Ik, and rrik are unit 
vectors; the "static" couplings between S and the {cJk} are described by the last two terms, and the dynamic effects 
are given by the terms in the curly brackets. The longitudinal static coupling tells us the change in energy of ak 
when S flips between ±zS; thus the vector wji Ik represents the difference between the effective field acting on ak 
before S flips and that field after S flips. If the {(Tfe} describe nuclei, this field will be almost entirely hyperfine in 
origin, unless the {ffk} are produced by truncating the levels of a higher-spin nucleus (i.e., a nuclear moment Ik with 
I Ik \> 1/2). In this latter case the total field acting on the {(Jk} will also include ,e.g., quadrupolar contributions; it 
then follows that the field uj^ ruk will be non-zero. This transverse field is the sum of the fields acting on ak before 

and after 5' flips. 

The dynamic terms in ([4.l| ) ori gina te as follows. The phase $ ^ t^S + J2k ^k is the sum of the "Kramers phase" -kS 



(see Introduction - notice that in (4.1 ) we have not yet included cos vriS* into the definition of Aq) and a renormalisation 
'Y^j. (pk of this phase, caused by the {(Jk}- The term ak describes the amplitude for ak to flip, under the influence of S, 
when S flips; both (pk and ak partly derive from the unitary "transfer matrix" Tk describing the effect of a tunneling 
S on ak'- 



Xk 



r>=e-*/''^^-*(^Mxr>^T? IxD; (4.4) 







(4.5) 



Here Hint{T) describes the microscopic interaction between S and ijk] \ X/T) environment state before S flips, 

and I Xfc™) after S flips. The signs ± refer to the path traced by S whilst flipping. The term in and further 
contributions to (pk and ak are necessary because environmental spins also have an effect on the tunneling action for 



S (in exact analogy with the nondiagonal terms discussed in Section III). 

The derivati on o f (4.1) has been described in detail elsewhere [0. Here we shall simply quote the values of the 
parameters in (4T) that arise if we start with the microscopic hyperfine Hamiltonian 



1 N 

HHyp = H^{S) + -Y,^kS-Ik 



k=l 



N 



k=l 



(4.6) 



i.e., our usual easy-axis/easy-plane model coupled to nuclei via hyperfine couplings Uk] we get 
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nS 

ak/V2-- 




(</>fc = 0) 

(if LOk < f^o) 



(4.7) 



4 



— X 

z 



(4.8) 



The assumption in (4.7), that the hypcrfine couphng Wfc <C i^o, is almost always true. In any case a.k^ £,k, and (f>k 
can be calculated as general functions of ujk/^o, for any initial microscopic Hamiltonian. It is useful to bear in mind 
the physical meaning of the parameter ak in this regime; ak is the amplitude for ffk to coflip with S. Note also that 
more general interactions than (|4.6|) will produce further renormalisations of Aq. 



Working with (4.1) (which is of course nothing but a low-energy effective Hamiltonian, in the usual spirit which 
also led to (2.4) and (3.1)) allows us to bring out the 4 essential physical mechanisms operating in the dynamics of 



S. These are 

(i) Even when Ho = 0, an internal field bias e = '^k=i ^\^k ^"-^^ '^'^ ^- Typically for nuclei we expect the hyperfine 
couplings to be tightly clustered around a principal value cjo, with lOq ^ Aq. There will however be a spread SuJk^ due 
not only to inter-nuclear couplings (principally via dipolar or Nakamura-Suhl interactions, with 5ujk ~ 10^ — 10^ Hz) 
but also other couplings such as transfer hyperfine couplings. One may then define a "density of states" W{e) for bias 
e (Fig. 3); if /i = N^/^6ujk/uJo <C 1, this consists of Gaussian peaks 



G^(e-WoA/V/2) 



exp{-(e-a;<,A/V/2)/r2}, 



(4.9) 



of width 



= M^o , (4.10) 

around e = uJoISN/2, where ZW ~ — Ni is the total nuclear polarisation, inside a Gaussian envelope of width 
N^^'^ujo] this envelope extends out to a maximum bias ±NuJo- 

W{e) = 2-^ Yl <^^^^'G,(e - c.„AA^/2) . (4.11) 

AN 

(Here we assume that N =even integer; otherwise AA^ = is impossible! The modifications required for N =odd 
are trivial and irrelevant to the ensuing discussion). In most cases, however, > 1, and the different "polarisation 
groups" (i.e., the different peaks) completely overlap, and we end up with a Gaussian envelope for W{e), i.e., 

W^(e)^— exp{-2eV(^o^)}; (m > 1) • (4-12) 

LUo 

where / = ^f^pKN . 

Now for tunneling to occur at all, the total bias (e -I- ^h) must be < Aq, otherwise the grain is simply trapped. 
The fraction of grains in an ensemble having bias < Aq is roughly A ~ /S.o/N5uJk if /i <C 1, and A ~ Ao/iV^/^Wo if 
/i > 1; the result for /i <C 1 arises because only a fraction ~ Ao/N^^^SiVk of that portion / = ^J^JtTN of grains having 
ZW = can fiip when the different polarisation groups do not overlap. 

Thus only a small fraction A of spins are not "degeneracy blocked" from tunneling (in the absence of nuclear 
dynamics, the effects of which are discussed below and in a Section 0). This de generacy blocking mechanism operates 
just as easily whether there is an external bias or not - in both cases only grains in the small window of bias can make 
transitions. 

(ii) Even in the absence of an external thermostat, the nuclear system can change the internal bias e. This occurs 
via dipolar interactions between the nuclei, at a rate ■ This pairwise flipping of nuclei conserves ZW but allows 
the nuclear bath to "wander in bias space" over the full range of bias energy associated with a particular polarization 
group, i.e., over an energy range ^ [ilOo — N^/'^SuJk- If the spread 5ujk is due entirely to dipolar interactions, so 
5uJk T2^, then the fluctuating bias covers the whole energy range in a time T2. This fluctuating bias e{t) can 
destroy coherence, but it can help magnetic relaxation, by helping the system find the small bias window. 
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(iii) The third mechanism, called "orthogonality blocking", arises because the transverse fields acting on the 

environmental spins (which are typically <C wjl) cause a mismatch between the initial and final nuclear wave-functions, 

in a way reminiscent of Anderson's catastrophe Defining sin2/3fc = w^/wj! (assuming j3k <C 1), this blocking 

effect can be parametrised by k, where 



N 



e-^ = \{cos(3k, (4.13) 



fc=i 



so that K « l/2X]fe/5fe- Orthogonality blocking inhibits both coherence and relaxation of 5; it spreads out the 
high-frequency (~ Aq) response of the grain to lower frequencies. 

(iv) A final interaction mech anism arises simply because the giant spin can flip the nuclear spins. This process is 



parametrised by the at in (4.£); for the whole spin bath one finds that on average a number A of spins will be flipped 
each time S flips, where A = l/2X]fcQ^fe fo^" <C 1. This causes phase decoherence; we have called this mechanism 
"topological decoherence" because it adds a random winding number to the effective action for S. It has a decoherence 
effect on tunneling, and thus inhibits magnetic relaxation. 

The formal discussion of all these mechanisms has appeared in our previous work. In the next 2 sub-sections we 
will deal first with the dynamics of S in the absence of spin diffusion between the nuclei, to show how the nuclei keep 
all but a tiny fraction of grains frozen; then we show how spin diffusion changes this picture. 

B. Tunneling Rate in a Bias 

We assume that neither the grain nor the spin bath is connected to an external thermostat except at time t < (in 
reality both are coupled to phonons and possibly electrons - see next section). We also neglect nuclear spin diffusion 
here, and estimate its role later. Thus relaxation of S, in a bias, can only occur if the bias fleld due to the {ak] allows 
states Ifl"; {f?™}), with the {ct™} in some polarisation state AA^™, to overlap in energy (within ~ Aq) with some states 
14; {f?fe*"}), where the polarisation state AN^™ ^ AN™ in general. Thus we expect some nuclear spins to flip, if S 
is to relax, although if /i is large, and the applied bias is small one can even have energy overlap between initial and 
flnal states having the same AN ^ (this does not, of course, mean that no spins are flipped during the transition, 
but only that none have to be flipped). 

The formal treatment of this problem is a generalisation of our treatment of the unbiased case. We write, for the 
time correlation function P{t) = {fz{t)Tz{0)), the form |l^ 

P(t; T;^h)^1+ [ de ^zm^^ ^ t^*'^*' ' + ~ " l] ' (^.14) 

This formula is crucial, and so we now spend a little time explaining it. The average J de is over initial bias, with 
both thermal and density of states weighting. PM{t,e) is a function deflned in Ref. jl^; it describes all those grains 
in an ensemble having bias e and for which, every time S flips, the polarisation state of the {crk} changes by 2M: 

PMit)= Tdye-y f2 [^Fy{iy)e'-<'''-^^Pl^,\t,e,^,y), (4.15) 
in which pfj describes a simple biased 2-level system (cf. Eq.(p7^)) 



,(0) _ -I ^lii^.v) 



= 1 - =5^Bin^(^M(^,y)0 , (4.16) 



El = e^ + Alj , (4.17) 



Am^v.v) - 2A„ I cos^Jm(2V(A- A')?/) | , (4.18) 
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with A = l/2X]fcC«fe; ^' = 1/2 Z]fe Q!fc('^fc)^ (so th^t ^ •^')- Eq.( 4.15 ) for PA/(t) can be understood as combining a 
"phase average" 

E / gFvMe^-(*-^) ; = e-^^--^ , (4.19) 

over the phase ip and winding number v), with an " orthogonahty average" 

dye-y (4.20) 



These averages are performed over the biased 2 -level correlation function P^^^ {t,t,ip,y), in which the tumichng 



ampUtude Am(</', y) depends on Af, Lp, and y via ( 4.18| ). 



Eq.( [4.14[ ) for P{t) may seem a httle odd, since it sums over apparently independent processes in which the polar- 
isation state of the nuclei changes by 2M each time S flips. Why should the polarisation state change by 2M each 
time? After all, when S flips, the number of nuclei which flip is random (with an average number A)! The answer to 
this question comes from energy conservation considerations. 

Consider a single grain that at t = has its nuclear environment in the state | xi,ei,zW) corresponding to the 
polarisation group ZW and the bias energy ei (see Fig. 4). Now 5 flips - what are the possible final states for the 
combined system, which have energies in resonance with ei? There are three possibilities: 

(1) None of the environmental spins flip; then the final state is | xii^/:"^^)) with e/ — — ei, since r^trl changes 
sign. 

(2) Some small number of environmental spins, say r <ti N, flip together with S, so that the net change in the 
nuclear spin polarisation is 2M (i.e., (r + M)/2 up and (r — M)/2 down spins are flipped). As a first approximation 
the final state energy will be « — ei + MlOq] assuming that the spread in nuclear frequencies is small, we may 
neglect for the moment the correction 5er ~ SuJkV^ (/zwo, Wo). We denote these states as | 2), and | 3) in Fig. 4. 

(3) An enormous number of nuclear spins flip together with S, with r as large as N] in this case we can not neglect 
the correction dcr^N ^ any more, because it is comparable or even larger (for /i > 1) than the typical difference 
Min[M}{'^£ ~ Mujo) ~ LOo- Then we can use 5er to "fine tune" a resonance ( which is impossible, in general, for 
the case (2)). Thus, we find in this case e/ = — ei + Mwo + Scr- The correction strongly depends, of course, on the 
particular set of nuclei flipped; a "fine tuned" situation is shown as state | 4) in Fig. 4. 

It seems at first, that the latter possibility is the best we can do for tunneling in resonance. Recall, however, that 
the probability that any given environmental spin flips during the transition is very small {a\ <C 1). Now, if we are 
going to tunnel between | 1) and | 4) then the amplitude of such a transition will be ~ Ao(a^)^, and for large TV will 
be even less then 2~'^ 5u)k\ Thus such a resonance is simply impossible to realise. 

In fact, only a small number of environmental spins may be flipped with a reasonable probability. Since the 
probability to flip none is just e^'^, we find the probability to flip exactly r spins to be pr ~ e~^{alYC^ ~ e^^X^ /rl, 
which peaks at r ^ A. The parameter A may be large in some systems , but still X ^ N. One immediately recognizes 
that the case (1) plays a role at small A, while the case (2) will dominate when A ^ 1 (formally, we could include (1) 
in (2) as a particular transition with r = 0). 

Now we make use of the inequality lUo 2> Aq to notice that among all possible transitions with different M, the 
energy mismatch between the initial and final states ei — e/ = 2ei — Mujo is either much larger then Aq for all M, 
or is close to the resonance ei — £/ ^ Aq for only one specific value M^j^, with all the other transitions being > ujo 
away in energy. Once we have a resonance the system will make transitions from to IJ- by changing the nuclear spin 
polarisation by exactly zL2M^^ each time to maintain this resonance. As for all th e oth er transitions with M =^ M^^ , 
they give only very small corrections, of order (Aq/wq)^ or less (cf. Eqs.(L7) and ( 4.16| )) to the main contribution to 



the dynamics of 5, coming from the M = M^^ term. Of course Eqs.(4.14) and ( 1.15| ) simply sum over all processes, 
and the resonance value of M is included automatically 

It is worth noting that nothing depend explicitly on the initial polarisation state. This point is further illustrated in 
Fig. 4. For given ei the initial polarisation state can be M, and if ei is close to the center of the Gaussian G^(e— Ma;o/2), 
then the resonance value is M^^ = M . If however ei is, say, the down-tail state of another polarisation group M + 2, 
then transition to the up-tail state oi —M + 2 has exactly the same energy mismatch and polarisation change 2M . 
That is why the only relevant statistical average is over the bias energy. Of course, if /i ^ 1, and different polarisation 



groups do not overlap, the bias energy is related to the specific group and Eq.(4.14) may be written as {S^h = 0) 



= 2^^ ^ ^N^''^'^ J deG,{e - Mojj2)e~P^PM{t. e + in- Mu^j2) . (4.21) 



15 



It is easy to understand now that external bias, which is indistinguishable from the internal one as far as the giant 
spin is concerned, will simply shift the resonance condition for given e to some other value M^^^^ , as is clearly seen 
from the diagram in Fig. 5, which is similar to Fig. 4, but now with S^h 7^ 0. 

We may further simplify Eq.(4.14) by changing variables from e + C-ff ~ Mujo/2 — s- e, and write it as 



M 



Z{(3) 



■[PM(i,e)-l] . 



(4.22) 



Then because the function e) — l] ~ (Ao/e)^ for large bias, and in most cases W{e) is a smooth function on a 

Ao scale (for ^lOq ^ Ao), we have 



P{t-T-^H)=Y. 



W(Mujol2 - Cff )e- 



-l3{Muia/2-i„) 



dtPM{t,e) , (/xWo>Ao) 



(4.23) 



(we assumed here that (3/S.o <^ 1). In the opposite limiting case fi = 0, when the Gaussian in (4.21) is a (5-function, 
we have 



Pit) = ^4(^E^i^^"^''^"''"°''^M(t,o) 



(4.24) 



M 



To summarize, we may understand (4.14) in a fairly simple way as including all effects of the spin bath on the 
dynamics of S. These come from the averages over phase (Eq.(gl9|), over the orthogonality mismatch between initial 
and final spin bath states (Eq.( 4.20| )), and over the internal bias e, acting on S, caused by the spin bath. It also comes 
from summing over all possible changes in the polarisation state of the bath when M flips. It is perhaps worth noting, 
for those who may be used to the theory of "oscillator bath" environments |17|, that the averages appearing here for 
the spin bath environment are very different in form, for the simple reason that in the case of the spin bath, most of 
the dynamics of the spins comes from their coupling to the macroscopic system itself {in our case, to S). By contrast 
in the oscillator bath models, the coupling to the macroscopic system is weak 0{N~^^'^), for each oscillator), and 
the dynamics of the individual oscillators is only weakly perturbed by the system. This is why, in the theory of the 
oscillator bath, one can first calculate the weakly perturbed oscillator motion as a function of the system coordinates, 
and then integrate out the oscillators by functional averaging. No such manoeuvre is possible for the spin bath (or 
indeed any other environment where the couplings are not weak), and it is not possible to map the spin bath onto an 
oscillator bath (see also Ref. [^). Any attempt to do so (by, e.g., writing a spectral function J{lo) ^ Nd{uj — tOo), 
for a set of "nuclear spin oscillators", as in ref. [^), cannot be expect ed to give meaningful results. 

From these remarks we can also see what has been left out of Eq.( 4.15 ) for PM{t); we have left out everything 
coming from the independent dynamics of the spin bath (independent, that is, from S). Formally this can be put in 
by incl uding these dyna mics in pj^"^ in ( 4.16| ). As we argue below this corresponds to allowing a time- dependent bias 
e{t) in ([4.14|) and (|4.15|) . 

We neglect the effect of coupling of the nuclei to phonons which is actually an incredibly small effect in most cases. 
In reality the nuclear spin dynamics comes almost entirely from nuclear spin diffusion, parametrised by the transverse 
relaxation time T2 , and caused by the dipolar interaction in ( |4.3| ) . As we shall see presently, if the spin bath is mainly 
composed of nuclei, this spin diffusion is crucial, since without it S cannot relax at all. Spin diffusion processes, in 
which AN remains unchanged but pairs of nuclei flip, cause e{t) to fluctuate in time because of the variation SiOk in 
the coupHng of each nuclei to S (thus, for a process |ti) ^lit) involving 2 nuclei, the total change in internal bias 
~ SiiJk', if N nuclei flip, the change in e is ~ SuJkN^/^). Now this fluctuation is fast, in the sense that for large values 
of AI the time it takes e{t) to change by an energy Am is usually much less than the time Ajj required for S to 
flip if it is in the coherence window (of energy width Am)', this will be demonstrated below. Consequently it is not 



necessary (or even useful) to go through the elaborate detour of recalculating in Eq. ( i.li ) , including a coupling 
to some effective " oscillator bath" , intended to model the effects of spin diffusion. We simply write instead that for 
this fast diffusion problem, where e{t) varies over an energy range fiujo, one has: 



P{t;T;^H)-l 



E 

M 



i^{N+M-Mh)/2 piM-MH)<^a/2 



2^Z{p) 



{PM{t,e{T)+S^„)),^r)-^ 



(4.25) 



where we introduced the notation = MhuJo/2 + S£^h with M^f=integer and | S(^h |< to define the shift in the 
polarisation change enforced by applying the external bias. Here (as before) 
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{PM{t,e{T) +6^h)) eir) = / dye 



E 



(4.26) 



(Pi?(i,e(T)+5eH,^,y))e(r) 



-t/TM(SiH) 



(4.27) 



(4.28) 



and A]\i{ip,y) is given by (4.18). In the usual case where > 1, or if the external bias is close to a multiple of Uo, 
the expression for the effective relaxation rate simplifies to 



' M 



(4.29) 



This concludes our formal discussion of the problem of quantum relaxation of S when it is coupled to a spin bath. 
We now use the results to find P{t,^H) for a few interesting cases. 

We start by ignoring spin diffusion, in order to demonstrate the way in which relaxation is blocked in its absence. 
Let us first consider a really pathological case, in which there is no degeneracy blocking at all, i.e., /i = 0. Then 
the function W{e) simply becomes a set of sharp lines, and all states in a grain ensemble may resonantly tunnel, if 
S£.H = 0. If S^H ^ 0, then no resonance is possible. Initially there will be fast relaxation, involving processes where M 
is not large. At longer times the higher-M processes take over - recall that for large M, the transition amplitude Am 
is very small, since Am ~ AoA*^/^/M!, which collapses when M ^ -s/A. To get some idea of the resulting relaxation, 
consider what happens if = 0, and let us ignore topological decoherence for simplicity (which makes no difference 
for the long-i asymptotics of the result). Then we have 



Pmit) 



dye-ypl^j\t,y) 



dy[l + cos{2AM{y)t)] 



(4.30) 



where Am(2/) — Aq JA/(2-\/Ay), assuming A' = (note if A' = 0, then the topological phase average ( 4.19 ) collapses to 
a (S-function S{ip — ^) and we get ( 4.3C ) anyway). The resultin g cur ves are shown for A = 5 in Fig. 6, for various PM{t)- 



It is easy to then find the behaviour of P{t) by substituting ( 4.30 ) into ( 4.24 ). The steepest-descent integral over M 
gives an accurate answer for long times. However we more quickly derive this long time behaviour from dimensional 
arguments. The sum over M contains some 



'2N polarisation groups inside the broad envelope function C^"*"*^'^^; 



thus the crossover from P{t) w 1 to the equilibrium P{t) w 1/2 occurs around a time t^ such that A 



'27V 



tr 



1, i.e. 



H^otc 



N f2N 
2-^Ha? 



For shorter times the relaxation is roughly logarithmic, viz. 



1 - P{t) 



2 In(Aot) 



{t < ic) , 



whereas for longer times one has 



P{t) - 1/2 ~ exp<^ - 3 



M ^'\~Aot) 

iVl/3 



[t > t,) 



(4.31) 



(4.32) 



(4.33) 



with logarithmic accuracy. The numerical evaluation of P{t) is shown in Fig. 7. 

However even the tiniest degeneracy blocking will upset these results - in fact if 5ujk\/N > Am, the result (4.3C) 
fails completely. Since the proverbial "thunderstorm on Jupiter" is enough to give a SuJk exceeding Am/^/N for the 
large values of M governing long time relaxation, we see that for relaxation (just as for coherence [0), it is the limit 
of strong degeneracy blocking that is experimentally meaningful. Let us again consider the zero bias case in this limit, 
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where > Aq/wo- We go back to (4.23), again assume A' = for simplicity, and, noting that the oscillatory sin 

d(0) 
M 



term in Pj^^ (Eq.(4.16)) gives an integrated contribution 



de 



dzJo{z) ; 



we find 



(4.34) 



(4.35) 



This result describes an ensemble of grains in which a fraction A of grains relaxes, leaving a fraction \ — A completely 
unrelaxed in the infinite-time limit; A is given by 



A^l-P{t^<:^)^^Y.— %(^^ dy\JM(2^) 



M 



(4.36) 



Since terms for M > ^/\ make almost no contribution to A due to the collapse of A^/, it can be approximated for 
large A by 



UJn 



(4.37) 



where / = ^^jixN is the number of states in the AA'' = polarisation group (in Ref. ||l^ we showed that individual 
terms in the sum over M decreased as \~^^^\ but the sum over all M gives the factor A^^^). 

Thus we see that apart from the small number of grains in near resonance, most grains will be frozen for eternity 
by the nuclear bias field, even in the absence of an external bias (again, ignoring nuclear spin diffusion). Adding an 
external bias makes essentially no difference to this (since this bias is physically indistinguishable from the internal 
bias). In fact the bias will have no effect at all on (4.35) until ~ the energy scale of variation of Wit). For 
/i < 1, this means that we will see oscillations in the decay rate as a function of with period ujo. For /i > 1, 
the decay rate and decay function A will change very little until ~ uJoVN; for ^ ijJoVN the relaxed fraction 
A — > since not even the nuclei can bring the system into resonance. We derive also a useful formula for the 
bias and temperature effects by noting that only small M <C \A/V contribute to the sum in Eq.( 4.23| ) and thus 
W{MlUo/2 - ^h) exp{-/3(Mw<,/2 - ^h)} ~ Wi-^n) exp{/3^j^)} giving 



WiO)Z{(3) 



(4.38) 



None of the results ( 4.30 )-( 4.38 ) is physical at very long time, because we have ignored the time vari ation o f e cau sed 
by nuclear spin diffusion. Let us now demonstrate that nuclear spin diffusion is fast, so that equations (4.26)-(4.28) are 
justified. Notice first that since dipolar pairwise flips occur at a rate for a given pair, and assuming short-ranged 
pair flips, we have roughly NT2^ nuclear flips per second in the sample. Then in a time A^/ we have N{T2A 



such flips, and in this time the bias will change by Se ~ StOo \/ NPiI^m ■ The condition for fast diffusion is 8t ^ Am 
(so the system has no time to tunnel), i.e., fast diffusion requires g 



^2 



(4.39) 



If the spread SuJo arises from the dipolar interactions between the same kind of nuclei only, then SuJo ^ (but 
usually SuJo > T^'^, because of other nuclei and fields), then this criterion becomes Aj^i T^^N^^^. Now in fact 
will be typical 10^ — 10® Hz (depending on the isotopic concentration of nuclear spins, etc.). For any mesoscopic 



spin Aq < 10® Hz (at least for most bare Hamiltonians), so that even if (4.39) is violated for M — (i.e., the zero 
polarisation change process), it is obeyed for large values of M. Since the bulk of the relaxation involves values of 
M ^ iV^/^, we see that (4.39) will always be the relevant condition. In this case Eqs.(4.25)-(4.29) follow immediately. 



Let us now therefore evaluate P(t) for the coupled grain/spin bath system, taking full account of the dipolar 
interactions between the spins, and assuming (4.39) to be satisfied; using (|05|)-(|427|), with A', this means we must 
evaluate, for ^ > Aq/ujo 
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Pit) -\=Wc'r 



-/3wo(M-Mh)/2 /-oo 



M 



dye ^ exp s — 2 



ri/2r, 



-tJit{2^J\y) 



(4.40) 



The integral is very easily evaluated; calling it /, we have 



dye ^ exp 



g Alt f\f_ 



M 



,M 



dye-ye{yM{t)~y) , 



where the step- function 6[yM{t) — y), with 



M2 



1/A/ 



arises because values of M range up to roughly ^/2N, and are thus typically large; thus 



M 



(n+m-Mh)/2 e' 



',3wo(Af-Mff )/2 



1 — exp 



l/Af 



(4.41) 



(4.42) 



(4.43) 



Noting again that if /i > 1, the dependence of P{t) on in (4.43) will be rather small until S.H ~ lOqVN , we begin 
by analysing this result for = 0. This zero bias case has of course been previously analyzed in our coherence papers 
[Q,^, but there we were only interested in the very short time behaviour of P{t), i.e., that part involving frequencies 
~ Ao, connected with possible coherent oscillations. Coherent behaviour can only arise in the M = polarisation 
sector of (4.30), which ignores spin diffusion. 

However in analysing the quantum relaxation properties of S, we need to sum over all ^/2N important terms in 
(4.43); the effect of the spin diffusion will be to unblock the long-time relaxation. In evaluating (4.40) or (4.43), a 
steepest-descent integration over M is possible, but just as with the case of no spin diffusion, dimensional analysis is 
sufficient; the crossover time from unrelaxed to almost relaxed behaviour occurs at a time tc given by 



In 



A2 



/27V In 



2N 



{AAA) 



(compare (4.31)), and the short-time relaxation is roughly logarithmic, as in (4.32), viz. 

— \n[2Klt/^^'^T, 



2nN 



In 



^ln(2A2t/^V2r, 



{t < te) , 



(4.45) 



The long-time behaviour is, analogously to (4.33), given by 

3 



P{t, ^H = 0)- 1/2 ~ exp ] - In^/^ (2Alt/7r'/'r^ 



{t > U) 



(4.46) 



with In In corrections in the exponent. Thus, amusingly, the results with spin diffusion included look just like those 
without spin diffusion, provided ^ — Q in the latter. The physical reason for this is simple - rapid spin diffusion 
basically eliminates the effects of degeneracy blocking, by allowing the bias e{t) to cover the whole energy range of 
each polarisation group. This is then the fundamental reason why spin diffusion "unlocks" S and allows it to relax. 

In Fig. 7 we showed a plot of P{t) against ln(2Aoi), which clearly brings out the 3 relaxation regimes. For very 
short times ^ A~^, there is a sudden relaxation involving only M ~ X^/"^ processes - the theory of this was given 
entirely in our coherence papers. For times t such that A~^ -^t-^tc, we have roughly logarithmic relaxation. Then 
for t '3> tc, we have the behaviour in ( [1.33 ) or [iAt). 

At this point it is important to realise that for all but microscopic spins, tc will be astronomically long. From (4.44) 
we have 



2A2 



2A^\ 



(4.47) 
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Consider now, as in Fig. 7, the case of a mesoscopic spin where, e.g., Aq ~ 1 MHz, = 100 MHz, N = 1000, and 
A = 10. We then find that tc ~ 10^° s, i.e., more than 10^° times the age of the universe! Thus for mesoscopic grains 
we always get logarithmic relaxation at long times, from nuclear spins, in the quantum regime. Readers familiar with 
the experimental measurements of magnetic relaxation in grains may well be quite surprised at this result, since it 
is almost a central dogma of magnetism that logarithmic relaxation must arise from a distribution of energy barriers 
pO| . Nevertheless we see that not only is this not so, but that the logarithmic relaxation has been derived here in a 
very general way, with essentially no approximations (but see below). 

On the other hand for microscopic spins {S ^ 10), one finds that tc can be short. For example, consider a situation 
where S = 20, = 10^° GHz, Ao ^ 100 MHz, F^lOO MHz, N = 10, and A = 0.1 (this would roughly describe a 
particle containing 10 Tb atoms). Then one finds tc ^ 1 s, and any direct relaxation experiment (typically conducted 
over time periods 1 ms < t < 10^ s) would be mostly in the long-time relaxation regime. 

These results are both surprising and interesting, in view of the very lo ng time (of ten logarithmic) relaxation 
observed in experiments. However we immediately point out that the results ( 4.45 )-( 4.47j) are by no means complete 



although they are correct to within In In corrections for a coupled grain/spin bath system, a realistic calculation 
must incorporate all three of the relaxation mechanisms we have discussed, in Sections ||, and this one. Since 
this changes the results yet again, we will not waste time here analysing the case of large bias {^h > i-^o\^ hi 



(4.25)-(4.28)) for the grain/spin bath system, but proceed directly to the general case. 



V. COMBINING THE MECHANISMS - PHYSICAL RESULTS 

We now come to the crux of the paper, which involves putting together the various mechanisms discussed in 
Sections ||-|^ to give a physically realistic picture of quantum relaxation in magnetic particles. As mentioned in the 
introduction, we do not attempt a comprehensive discussion of all possible cases. Such a discussion would involve 
consideration of a large class of "giant spin" bare Hamiltonians, with widely varying behaviour, and with each being 
considered throughout the whole range of field strength, field orientation, and temperature. This would only be the 
beginning - we would then have to go on to examine the large variety of couplings to nuclear spins (with hyperfine 
couplings ranging over 3 orders of magnitude, plus quadrupolar couplings), the variety of different magnetoacoustic 
couplings, etc., etc. 

We adopt the tactic here of showing how the mechanisms combine, and giving some details for one model. In other 
papers, in preparation, we analyse some specific experimental systems in much greater detail. We emphasize once 
again here that, in our view, magnetic relaxation in most experiments involves the dipolar interactions between grains 
in an essential way. 

We begin by showing how phonons and nuclear spins combine to give a rather surprising form for the relaxation 
at low T. This analysis is appropriate to the case where both grain and substrate arre insulating. We then go on 
to include electrons - in this case phonons become irrelevant, and only the combination of electrons and nuclei are 
important. 

Before discussing the details, let us first state the method we shall use. We have already seen how the dynamics of 



S is given in the presence of a spin bath (Eqs.( 4.25 ) or ( 4.4C )). When we couple in a bath of oscillators which allow 



transitions between states | S'i;xi{CTfc}) and | S2',X2{'^k}) of the combined grain/spin bath system, in which there is 
no restriction on the difference 2M between the polarisations of the spin bath states xi) or X2), or on their biases ei 
and 62 (apart from those imposed by energy conservation), the form of P{t) must change. In this paper we will not 
attempt to give a general expression covering all mechanisms simultaneously; this is more than we need. Instead we 
will use 2 expressions, which apply in the 2 limiting cases of interest. These are 

( a) Fo r short times, when the spin-mediate relaxation dominates; then we ignore the oscillator bath, and use 



Eq.([4.14D 



(b) For longer times the oscillator bath-mediated relaxation will take over. In this regime, each time S fiips, some 
of the bias energy is taken up by an oscillator mode, and the rest by some number r of fiipped spins in the spin bath. 
There is no restriction on either r or the polarisation change 2M, for this process to work. The complete calculation 
of P{t) here is very complicated, but we notice that as soon as we reach times where the typical oscillator-mediated 
transition rate, at a typical bias, is faster than the spin bath-mediated transition rate (for those grains in an ensemble 
which have not yet relaxed), then a reasonable approximation for P{t) will be 

P{t; T,^h)^J de^^^^^e-^^^-^"^ |p(^') (T, .) + [l - p(^^) (T, e)] e'^'^^^^^ | , (5.1) 
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where r ^{e,T) is the oscillator-mediated relaxation rate, and P^'''i\T,e) = e '^^'^ /{2cosh{e/T)) is the equilibrium 
population of the state Si in a given bias (in our previous discussion of the nuclear spin effects we assumed that this 
bias was actually much less than T). We note that this expression is incoherent, because it applies to the majority of 
grains in a large bias, which cannot relax via the spin bath as discussed in the previous Section. The relaxation rate 
in this incoherent case can be calculated as a second order perturbation theory expression in the tunneling amplitude 
Aq. Thus the only effect of the spin bath which is left is the distribution over the bias. 

The formal proof of the above statement is easy to follow for the case of pure othogonality blocking. In a typical 
bias e ~ WoiV^/^ 3> Aq the second order perturbation theory expression for the tunneling rate is given by 

r-i(6) = 2nJ2pr I U.f \' E z^"' I ^^"'^ I' + - E, + - Mf)/2) , (5.2) 

if "7 

where the sums are over the initial and final states of the spin bath (i, /) and oscillator bath (a, 7) with the equilibrium 
density matrices p'^'^ . Here describes the interaction with the phonons (we do not even need its explicit form to 

prove the point), and the orthogonality rotation operator U is given by 

N 

U =Y[ e-'^^*" . (5.3) 
fe=i 

The crucial point here is that in a large bias one may drop the energy LUo{Mi — Mf)/2 transferred to the spin bath 
from the 5-function argument in (|]^), because only transitions for which Mi — Mf ^ iV^/^ contribute to the answer. 
After that we have the sum over the complete set of states /, so that UifUj- = 1 because the of the unitarity of U, 
and the final answer is that for the oscillator bath alone. Including topological decoherence changes the result in a 
minor way. Indeed, with nonzero ak, we have 

N 

fc=l 

and the result for the relaxation rate in the oscillator bath is renormalized by the factor 

£> = (f/;/^) = 2 + 2cos(2$)Fa(1) , (5.5) 

Unless A and the topological phase $ is a multiple of vr, this factor is of order ~ 1, and clearly has no essential 
effect on the results. Thus, apart from the renormalisation factor D, the relaxation of a single grain in a nuclear bias 



field e is that already given in sections ^ and III, and all that remains for a grain ensemble is to average over the bias 
field. 



A. Spin Bath Plus Phonons 

The case of insulating grains and an insulating substrate is the one where there is most obvious competition and 
interplay, between the nuclear spin and phonon relaxation mechanisms. They are both slow. We have already seen 
(Fig. 6, and Eq.( 1.45| )) what the short-time relaxation will look like. Now consider the implications of ( pTl| ) for the long- 



time relaxation. We worked out r^^(e) for phonon relaxation in the absence of nuclear spins; for a longitudinal bias 
this was given by ( |3.13| ), and will be roughly the same in the presence of the spin bath (i.e., with the renormalisation 
factor (|5.5|) above). Thus we write r~^(e, T) in the form 



e 



3 



r-i(e,T) = T-i ( — ) coth(e/2fcBT) , (5.6) 



where ^ is a "typical" normalising relaxation rate, defined as 



3 



r-^^r-\e^E,MT^E„)^S^l^4^(^] , (5.7) 



and Eo is the width of the Gaussian peak in W{t), i.e.. 
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(5.8) 



Notice that for a fixed bias e = Eo, the actual relaxation rate r ^{Eo^ T) still depends on temperature as r ^(e, T) = 
r^^ coth(£^o/2fcBT). At low T this scales down to the constant value t^^ , whereas at high T it increases linearly, i.e., 
T-\Eo,T)^2T-^kBT/Eo. 

It is useful to get some preliminary idea of typical time scales for Tq. Let us consider 3 examples, viz., (i) a Tb Oxide 
grain, with S ^ 1000, containing N ^ 1000 nuclei; (ii) a particle of Er As (such particles are apparently insulating 
containing N - 1000 nuclei, with S ^ 1000; and (iii) Ni O grains with S ^ 1000, but iV ~ 10 only (coming from 
the 1% of Ni'"^^ nuclear spins). Let us also assume that Qd ~ 100 K and Aq 1 MHz for each example. However 
Eo varies dramatically - for the Tb O grain, Wo ~ 5 GHz, and Eo ^ ^ K\ for the Er As grain, Wo ~ 1 GHz, and 
Eo 0.8 K; whilst for the Ni O grain, ujo ~ 28 MHz, and Eo ^ 2.3 x 10"'^ K. We thus get a wide range of time 
scales; for the Tb O grains one has To ^ 1.6 x 10^^ s, for the Er As grains, one has To ^ 5 s, whilst for the Ni O grains, 
one has Tq 4.5 years! 

We now consider the general behaviour of P{t,T; ^h)- We start by noticing that the weighting function for the 
initial states, of bias e, contributing to ( ^.l|) , will be 



wie,T;^H) = W{e~CH)- 



-/3(e-?ff) 



Z{f3) 



1 



2ttE, 



■ exp 



Eo 



Eo • 



(5.9) 



Ignoring for the moment, we see that in zero external field this distribution peaks at negative bias e —E'^/ksT. 
This is natural; as we lower T, more and more of the nuclei align with S, thereby lowering the energy of the combined 
system in the initial state. This process continues until ^ N nuclei are aligned with S, i.e., when e ~ —NuJo/2, at a 
temperature fcsT ^ uJo/2. We have then reached the bottom edge of the distribution W{e), way out in the wings. If 
7^ 0, the distribution peaks at e ^ {S^h — Eo/kBT). 

Thus we write 

Pit,T;^H)^ J de^(e,T;eH)|p''=^HT,e)+ [l-P(^«)(r,e)]e-*/^(^'^)| 



dx 



/27r 



„a:/Tg-(t/To) coth(2:/2T) 



2cosh(x/r) 



(5.10) 



where we have normalised all energies by Eo- 



X = e/ Eo 

= in/Eo 
T = keT/Eo 



(5.11) 



Let us first look at the high-T limit of ( 5.10| ), i.e., for fc^T ^ Eo when T ^ 1. Then ( 5.10 ) is easily evaluated; one 
finds 



Pit,T;^H)-l/2 



1 1 

2(l+4Ti/To)l/2 

1 1 

2(l + 4Tt/T„)i/2 



exp 



2Tt 



{keT > Eo) 



(5.12) 
(5.13) 



(We recall that this result is correct only when T flo - otherwise higher levels of the grain will come into play.) 

This surprising power-law behaviour is actually easily understood. The initial fast relaxation (at t ^ To) comes 
from grains with large bias; the slower relaxation (t ^ To) comes from the grains with smaller bias. The typical 
relaxation occurs at a rate TTo^, i.e., it is faster for higher temperatures, in accordance with the remarks just after 
Eg. ( ^.81 ). We may think of the power law as a "grain-ensemble" sum of a lot of different exponential decays, or as a 
funny kind of "stretched exponential". Note however that the decay is not logarithmic in time. 



Now consider the low-T limit, with T <C 1; this is a little more complicated. From (5.10) we now have 
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dx 



,-i(x-C„ + l/T) = 



^ COt\\(x /2T)t/Ta 



'2tt 



2cosh(x/r) 



27r 2cosh(a;/T) 



coth(2;/2T)t/To _ -j^ 



(5.14) 



Let us consider first the behaviour of (5.14) when = 0. Using the fact that T <C 1 one may further simphfy this 
expression to 



T 



dx 



cosh(a:) 



This clearly defines a temperature-dependent relaxation rate 



-x^ cotli(x/2) tT /t 



e 



D 



(5.15) 



(5.16) 



which goes to zero at low T. We note however that the validity of ( 5.15 ) requ ires t hat the calculation is not affected 
by the edges of W{t), which means that at a temperature ksT ~ Eq.( 5.1(: ) crosses over to a constant value, 

i.e., 



\2eL 



(5.17) 



but the magnitude of the relaxing component at this temperature is already extremely small; in fact 1 — P{t — > 
oo) ~ 7V~^/^e^^. It is again useful to consider what is the the maximum relaxation time at low T. For the 
examples previously mentioned, we have (i) for the Tb O grain, Teff{T ^ 0) ^ 10 minutes; (ii) for the Er As grain, 
TeffiT ^ 0) ~ 2 days; and (iii) for the Ni O grain, Te//(T 0) 4 x 10^ years. 
Defining the exponentially small total amplitude of the relaxing component as 



A 

we find the long time asymptotics of P{t) as 

P{t,T;iH^O)- 



2V2 



l-A + A 



2reff{T) 



TTt 



1/2 



(5.18) 



(5.19) 



The reason for the decrease in the relaxation rate at low T is simply that as we lower the temperature, the typical 
bias energy decreases linearly with T, because only a small fraction of grains in the Gaussian tail which have very 
small bias e ~ fc^T contribute to the relaxation. The majority of the grains are actually "trapped" by the hyperfine 
interaction in a negative bias energy e ^ E'^/ksT ^ ksT. Thus at this low temperature even phonons will not help 
to liberate the giant spin. 

Consider now Eq.(5.14) with nonzero external bias. For <c 1/T the previous answer ( 5.19| ) hardly changes 
except that the amplitude of the relaxing component is given now by 



2V2 



i^H « 1/r) . 



(5.20) 



We observe that already for rather small bias ^j:^ ^ T <C 1 we have an exponential dependence on . The negative 
bias will simply further suppress an already exponentially small fraction of relaxing grains and we shall not go into 
more details here. The case of positive is much more intriguing and surprising. After an exponential increase in 
the amplitude up to A 1 for ^ 1/7, the answer changes drastically. For > 1/T one finds the magnetisation 
function to be 



dx 



(5.21) 



which has different behaviour depending on whether we look at short times, where 
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P(t,r;C«) « e-*/-=/H«-^) , (Vre// < (Ch - 1/^)') , (5-22) 

or long times, where 



r:h{^H,T) = (Ch - ; (5.23) 



P(i,T;e.) ^ ^^E^e^il.-yTfn {I^f' , (,/,^^^ > (J^ _ i/T)^) . (5.24) 
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Note that Eqs.(5.21)-(5.24) are derived in the approximation that — l/T ^ 1. As before, the bias is restricted to be 
much less than fio. Finally, at t > Tq, this behaviour changes yet again to P{t) w exp{ — (^^ — l/T)^/2} (To/4tr)^/^ 



We note the anomalous temperature and bias dependence of the effective relaxation rate ( ^.23 ) near the crossover 



From this analysis we conclude that unless the external bias is larger than Eq (or E^/ksT at low temperature) 
the ensemble averaged results for the magnetisation relaxation are not described by the naive theory of a two level 
system coupled to a bath of oscillators, and spin bath interactions essentially modify the answer. It is worth noting 
three effects which have to be kept in mind when considering the evolution of the experimental data in external bias: 

(i) Since the dependence on starts when > Eo or > E'^/ksT, one cannot derive the value of the tmmcling 
amplitude from this dependence, as it would be in the case of isolated two-level system; 

(ii) At low temperature the giant spin is "trapped" in its initial state by the large and negative bias produced by 
the spin bath. 

(iii) Only large and positive bias can liberate the giant spin and to allow complete magnetisation relaxation. The 
temperature and bias dependence of the relaxation rate is anomalous near the crossover region ~ 1 /T. Of course, 
when 3> 1/T we recover back the pure case of oscillator bath relaxation. 

At this point we would like to comment on our starting assumption that the initial state of the spin bath is 
equilibrated with the giant spin direction S — Si. Experimentally this could be arranged by applying a very strong 
negative bias S.H during a time period much longer than the longitudinal NMR relaxation time Ti , and switching it 
off at t = 0. In some systems however this procedure may not work because of an astronomically long Ti (in fact, 
low temperature longitudinal NMR is still something of an unsolved mystery; spin-lattice relaxation times at mK 
temperatures ought to be many years in insulating crystals, unless some gaplcss magnetic excitations are involved, 
e.g., on the sample surface or crystal defects). It may be useful then to introduce two different temperatures - one 
for the oscillator bath Tph and the other for the spin bath Ts (which may even be negative, i.e., Ts < 0!). Our basic 
equation ( |3.l[ ) is still valid in this more general case, but now the distribution over the initial bias is defined by the 
spin bath temperature, whereas the giant spin evolution toward equilibrium is governed by the crystal temperature, 
i.e., 

P(t,^H) = J &^^^|^^e-(-««)/^-|p(^'^)(rp,,6) + [1 - p(^^)(rp,,6)]e-*/^(^-^^")| . (5.25) 

One may proceed with the analysis of this expression as before. At this point we feel that considering more cases in 
this paper will not add much to the physical picture. Depending on the particular experimental system and sample 



preparation, the necessary formulae can be easily derived from (5.25) 



Summarising what we have found for insulating systems, we see that at short times, the relaxation proceeds entirely 



via the nuclear spin bath, and is logarithmic in time (Eq.(4.45)). At longer times phonons take over, and we get power 
law decay in time; at low T this goes as {T^f f{T)/t)^^^ , with r^.^ ~ T"^ (Eqs.( 5.16 ) and ( 5. IS )), but only a small 



fraction of grains relax unless the bias ^ Eo/ksT or larger. This bias is necessary to counteract the nuclear 



bias field. At higher temperatures we still get power-law relaxation (Eq.( 5.13 )). The crossover between nuclear 
spin- mediated and phonon- mediated relaxation can be understood by matching Eq.(4.45) with the relevant phonon 
expression. 

B. Spin Bath Plus Electrons 

We now turn to the case where either the grain or the substrate is conducting. We use the giant Kondo model of 
section ^ to describe interaction of S with electrons, with a dimensionless coupling a. If both grain and substrate 
are conducting, then a — ab ^ g^S"^^^ (Eq.( |2.7| )); if only the substrate conducts, then a = c^s ^ g^S^^^ (Eq.(2.13)); 
and if only the grain conducts, a ~ ase'~^^''^'^^ (Eq.(2.21)). 
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The formal analysis is almost identical to that just used for phonons - we start again from (|]^), now using the 
general result (2.8) for the electronic relaxation rate. The high temperature limit is most transparent because Te(T) 
in Eq.(2.9) is independent of the bias. Thus we find a pure exponential relaxation 

P{t) = 1/2(1 + e-*/^=(^)) ; {aT, T . (5.26) 

at high temperature. 

For very small grains or conducting grains on an insulating substrate at T <C we have a peculiar relaxation 
regime when the internal bias is much less than temperature, but the electronic damping rate = 2TTaT is already 
small, that is for a <C 1 there is the temperature range where F ^ E'o <C T. Now the electronic relaxation rate is 
inversely proportional to the bias energy 

A2 

r-i (T, e) = STraT-f , {2TrT < e < T) . (5.27) 
Substituting this expression to ( ^.1| ) we find the time correlation function as 

P{t,T;^H) + -^e-^(-«-)%-*/(-^-(^'^°) , (5.28) 
which after easy evaluation yields 



P(t) = 1/2(1 + e-VW^^(W) ; (^^ = 0) (5.29) 
« 1/2(1 + e-*/^=(^'««)) ; (Cff » Eo) . (5.30) 

The temperature dependence of the effective relaxation rate is Ta{T), and again in the small external bias we have 
an unusual decay law - it is neither simple exponential nor power-law. 

Let us deal now with the low-temperature behaviour, ksT/Eo <C 1. Using (2.1C) we then get 

F(t,T;^H)«l+ / ^e-i(^-«-+i/^)^-4f-=-(e-/(^)*/-//-l) , (5.31) 
J V27r 2cosh(a;/T) V / 

f{x) = cosh(x) I F[a + ix/n] p 
(compare ( ^.14 )). If the bias < 1/T, this defines a temperature dependent relaxation rate 



'2ttT 



2«-l 



(5.32) 



The answers crucially depend on the parameter a. We consider three limiting cases: 

(i) For a <C 1 we approximate the function f{x) by f{x) « a;coth(x)(a^ -I- {x/n)'^)^^ and obtain 

Trf 

P{t,T;^H)^l-A ; (i«aVe//), (5.33) 



TTt 



1/2 



P{t,T;^H)^l-2A[ ; (a'T^f j < t < r^f f) , (5.34) 

/ f , 

P{t,T-in)~\- A + a[ e-2V*/-.// ; (r,// » t) . (5.35) 

(ii) For a = 1/2 we have /(x) — tt, and, as in the high temperature limit, the relaxation is given by the simple 
exponential law 

P{t) = 1 - A(l - e-^'l^^s!) ; (5.36) 

(iii) Finally we consider the case of large a, which has a power-law asymptotic at long times 

On /T c t \ 1/(2q-1) 

P(,,r;M»l-.4 + .42r[^](iai) ,„7, 
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In all these expressions A{T,^h) is given by (5.20) 



The interpretation of these results also depends on the value of a. If both grain and substrate are conducting, then 



the relaxation time in (5.32) will be greater than the age of the universe unless S < 100 (depending on the value of 
g^). This shows the astonishing power of the electron bath to "freeze" the dynamics of S, unless it is quite microscopic 
in size. We shall not analyze here the case of very small S, since once S ~ 0(10 — 100), the detailed structure of W{e) 
will no longer be Gaussian - there will be lots of fine structure, depending on the particular systen involved. Such 
studies are best done on a case-by-case basis. 

If only the substrate is conducting, then a — as will be small until S ^ 10'^. For larger S then the giant spin is 



again frozen. For S < 10 the results (5.33)-(5.36) can be applied if S,ff ^ 1/T. Again, as approaches 1/T, there 



will be an amusing crossover , with a relaxation rate given for finite S,^ > ^/T by 

t;/^. (T, ^h) = (T, e = ^H- El/T) , (5.38) 
Finally, if only the grain is conducting, a can be exponantially small, if the spacing between the internal electron 



levels is large enough - this typically requires that Ae^ > Vlo, so that S < 10^ (cf. Eq.(2.23)). In this case, unhke the 
cases where the substrate is conducting, the short-time relaxation will be again dominated by nuclear spin-mediated 
transitions, and the weak electron-mediated transitions only enter at later times in possible competition with the 
phonon- mediated transitions; the typical electron-mediated rate becomes 

r-)iT) = inas^e-^^/^ , (5.39) 

Note that for low enough temperature, r^ff in ( |3.39| ) becomes so long that electrons in the grain become irrelevant, 
and we return to the calculations for the nuclear spin/phonon- mediated problem, i.e., the grain behaves as an insulator. 

Thus to summarize, we see that depending on whether the grain and/or substrate are conducting, we get behaviour 
ranging from grains frozen for all eternity (if a 3> 1), to grains behaving essentially as insulators, with short-time 
relaxation controlled by the nuclei (when a — > 0). In the former case we see an extraordinary alliance between the 
electrons (with their strong dissipative suppression of tunneling) and nuclear spins (which trap in a negative bias 
e ~ —E'^/T) to block any motion of S even at low T, where traditionally one expects tunneling. In fact, as T ^ 0, 
our calculations show that S is blocked for eternity when S takes microscopic values (e.g., S ~ 0(10)), because the 
negative bias becomes very large. This is a quite astonishing demonstration of the control that nuclear spins exert over 
the dynamics of S - the spin is only unblocked when S ~ 0(1), when a <C 1 again (and the whole WKB framework 
breaks down!). 

VI. SUMMARY & CONCLUSIONS: PRACTICAL IMPLICATIONS 

Let us first summarize our main results. At temperatures T such that fc^T ^ fJo, we may employ a 2-level model 
with bias to describe the grain. Coupling to electrons and/or phonons in this model can be treated by the usual 
oscillator bath methods. However the crucial coupling turns out to be to the nuclear spin bath, to which oscillator 
bath models do not apply. The main effect of this coupling is to spread each of the 2 levels for S into a Gaussian 
"band" of half- width Eg ^ N^^'^ujo, and total width Nuq, where lOq is the hyperfine coupling and N the number of 
nuclei in the grain. If TV and oJq are small, Eq may only be a few mK (as in Ni or Fe grains). If N and lOq are large 
(as in mesoscopic rare earth grains), Eo may be hundreds of Kelvin or more (so that i?o ^ ^lo)- At low temperatures 
the nuclei begin to line up with the grain vector S*, and then S finds itself in a negative internal nuclear bias field, 
with a mean bias e of roughly e ~ —E'^/T. Then at short times the only way that the grain can relax is by taking 
energy from the nuclear system. The inter-nuclear dipolar interactions allow this to happen, by causing the bias e to 
become time-dependent - the result is a slow relaxation, which for an ensemble of grains gives a fraction of relaxed 



grains going roughly logarithmically in time (section IV) 

1 - P{t) ^ 









2ttN , 
in 


1 

e\/A 


In (2A2i/7ri/2r^) 





(6.1) 



This logarithmic relaxation, shown in Fig. 7, has nothing to do with a distribution of sizes of the grains (we assume 
all grains have the same parameters) . 
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At longer times the grains can relax very slowly via either electron- or phonon-mediated transitions, again in the 
time- varying nuclear field. If the grain is insulating (no electrons) , we get power-law relaxation; a fraction A of grains 
can relax, where 

A « exp ( - 1/2 f 1^ - ^ V I ■ (6-2) 



2Eo I \Eo T 

in a bias field ^h, and these grains relax as A{Tef f /ty^^ where 



3 



r-/,(T) ^ 5^aJ ^ ) f ^ ) , (6.3) 



if < E'^/T, and P{t) - 1/2(1 + e^*/^-" ) with 



?2 



once > /T . Thus very few grains can relax at until the external bias S^h compensates the internal bias e and 
"untraps" S. 

If there are electrons around, the results depend on whether the substrate is conducting or not. If it is, then all 
but microscopic spins (S* <C 100) will be trapped in states of frozen magnetisation for astronomical times, once the 
temperature goes below VIq. This remarkable result is a combination of "degeneracy blocking" caused by the nuclear 
bias field, and the very strong electronic dissipation. However if only the grain is conducting, and the substrate is 
insulating, then for S <C 10^, the electronic relaxation becomes negligible, and we go back to the fully insulating case; 
the full results for all values of electronic coupling appear in section 0. B. 

It is clear that these results have implications both for future magnetic device design, and for the thousands of 
experiments that have been done on magnetic grain relaxation over the years. As far as devices are concerned, 
perhaps the most interesting result is the freezing of 5* when both grain and substrate are conducting unless S <ti 100. 
This may have far-reaching implications for future computers and information storage, since it implies that if we are 
prepared to go to low temperatures, one should be able to use even magnetic molecules as permanent memory storage 
elements, perhaps of the "M.R.A.M." type (cf. ref. |^l|). Even more tempting is the possibility that we may be able 
to manipulate S in such molecules indirectly, by controlling the nuclear spin polarisation, since it is the nuclear bias 
field that controls the dynamics of S. 

As far as experiments on grain relaxation are concerned, we believe that it is quite urgent that experimental tests 
of our results be done at low T. Unfortunately most of the many relaxation measurements that have already been 
done are at higher temperatures, where relaxation is dominated by thermal activation. What we have found here 
is that the low-T behaviour does not look like the conventional tunneling picture at all - there is no crossover to a 
temperature - independent "tunneling relaxation" at temperatures fesT <C ^o- The low-T behaviour of magnetic 
particles is not governed by the external bias at all, unless it is very large; instead it is governed by the random 
distribution of internal nuclear bias fields. This renders previous theories of the tunneling of magnetic particles 
somewhat meaningless, regardless of the external field. Thus we believe that it would be useful to look again at those 
experiments that have been done at low T on ensembles of relaxing grains [p[-p| j2^j30[ | , particularly in an attempt to 
understand the "plateau" that often appears at low T in the relaxation rate (the "magnetic viscosity plateau"). 

However an important cautionary note is necessary. In virtually all of these experiments, inter-grain interactions 
(mediated by the grain dipolar fields) are very important. We believe this to be the reason for the common occurence 
of "avalanche" magnetisation reversal in grain ensembles (see, e.g., ref. |2^ ); these avalanches are prima facie evidence 
for the importance of inter-grain interactions (such avalanches also occur in multi-domain magnets, because of dipolar 



interactions between the domains and walls |7|j32[|). Thus to interpret most experiments (including the remarkable 



recent results on Mni2 0i2 molecules 22 1), we must incorporate these interactions - this we do elsewhere [ p3[ . 

However some experiments may not suffer this problem. The recent results of Coppinger et al. |3^], on ErAs grains, 
with S ^ 10'^, may be a case in point. It is not yet clear what is the effective coupling to the electron bath for this 
system, but we believe that the theory in the present paper should be applicable to this case. 
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APPENDIX A: KONDO-TYPE RENORMALIZATIONS IN THE MULTI-ORBITAL SCATTERING 



In truncating the initial Hamiltonian (2^) to the 2-level system for the giant spin one generates both the diagonal, 
i.e., proportional to TzC^z term, and nondiagonal terms ~ T±az^. The latter are usually not present in the spin-boson 
Hamiltonian because they are subdominant to the diagonal ones. In general we may write 



Ck'l' 



(Al) 



kk'l 



kk'U' 



Here I and I' are the orbital channels for electron scattering off the grain {Imax ~ kpRo ^ 1), and we write the 
Hamiltonian to be diagonal in I for f^fj^ coupling. We observe that initially (before scaling) the transverse coupling 
J-^ is exponentially small as compared with J^. In fact J-^/J^ Ao/ilo ~ (tfl-IJ-) because transverse matrix elements 
flip the giant spin as well and are therefore proportional to the overlap between the Iff) and |4) states. 

Formally (Al) is a strongly anisotropic Kondo model with multi-channel scattering in a magnetic field Aq applied 
along the weak-scattering direction. Since the J"*" are small we may ignore the renormalization of the and study 
the renormalization group equation for the first. Introducing dimensionless couplings gl[, = J^,N{0) and the same 
for of, we find 



9w 



9w 



Si + Si' 



(A2) 



1 

2^ 



(A3) 



where Si is the scattering phase shift due to . The crucial point is that for a large grain, when many scattering 
channels are equally coupled to the giant spin, the Anderson orthogonality term as will always dominate in the 
scaling of g'^. Thus, quite generally, we never face the strong coupling regime when the g-^ are of order g^ and large. 

Naively one would ignore these transverse couplings completely, but we have to remember that the giant spin 
dynamics is governed by the transverse magnetic field Ao which is also subject to renormalization due to Anderson 
orthogonality. The relevant renormalization group variable is fugacity Ao/T, and 



A. 



A^(l-a^ 



(A4) 



(the solution of this equation is explicitly used in our equations (2.8)-( pTTc| )). Thus the question of whether we may 
neglect conduction electron spin flips or not depends on the parameter 



e = Max{^^} 



(A5) 



If is less than unity we are back to (pT^)-( 2T0| ), and electron spin flips are irrelevant. If we find > 1 then the 
correct expression for the giant spin relaxation time will be 



27T I 4 1 2 7V(0) 



(27rT, e)r 



s,, 



(A6) 



Note that for 6* > 1 we definitely have as '-^ S"^/^ 3> 1, so the S dynamics are definitely incoherent, and the above 
refinement is hardly different from (2.8)-(2.10) (it is equivalent to going from 2a.s — 1 to 2as — 0). 



FIGURE CAPTIONS 
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Figure 1 The Atomic Weight (in Dahons) of a typical computer memory element, as a function of time. Note that 
the position of the "Quantum Threshold" depends on temperature. As we see in this article, it also depends strongly 
on what the element is made of. 

Figure 2 Schematic diagram of the energy scales in the giant spin, biased by a longitudital field Hz- The bias 
e = 'feSHz, and Qo is the gap to the next set of excited states. In reality, for a giant spin, the top of the barrier will 
be much higher, at energy ~ (Silo) above the lowest states. 

Figure 3 The distribution function W{e) for spin bath levels, around each giant spin level. Typically the parameter 
/i > 1 except for very small numbers of nuclear spins (see text). 

Figure 4 Some possible transitions when S flips, for a grain in zero external field. The initial state of the combined 
giant spin/nuclear spin system has energy ei; in the diagram this state | 1) has either initial nuclear polarisation 
AN = M 01 AN = M + 2. Transitions to states | 2) or | 3) are accomplished by changing the nuclear polarisation 
by 2M (so the final polarisation is either —M or —M + 2). The transition to state | 4) is made by flipping a very 
large number ~ of nuclear spins, whilst still changing the polarisation by 2M only. This allows us to "fine tune" a 
resonance with state | 1), as described in the text. Only a few polarisation groups are shown in the Figure; the insert 
shows how these fit into the distribution W{e) (shown for both initial and final states). 

Figure 5 The same set of transitions as described in Fig. 4, but now for a grain in an external bias £,h, acting on 
S. The change in polarisation is Mi + M2- The insert shows the initial and final distributions W{e zt £^h), displaced 
from each other by 2(^h- 

Figure 6 The time dependence of different contributions Pnit) to the time correlation function Pit) for an 
ensemble of grains interacting only with nuclear spins, ignoring the effects of nuclear spin diffusion. Wc assume that 
the parameter A = 5, i.e., roughly 5 nuclei out of A = 1000 are flipped each time S flips. 

Figure 7 The time correlation function Pit) for an ensemble of grains, now plotted against Int; all contributions 
Pmit) are included. This includes only nuclear spin bath- mediated transitions. The figure describes either (a) P{t) 
neglecting nuclear spin diffusion, and having zero degeneracy blocking (^ = 0), or (b), including both spin diffusion 
and finite degeneracy blocking. 
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